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Abstract 

In this paper we use solar neutrino data to derive stringent bounds on Majorana 
neutrino transition moments (TMs). Should such be present, they would contribute 
to the neutrino-electron scattering cross section and hence alter the signal observed 
in Super-Kamiokande. Motivated by the growing robustness of the LMA-MSW 
solution of the solar neutrino problem indicated by recent data, and also by the 
prospects of its possible confirmation at KamLAND, we assume the validity of this 
solution, and we constrain neutrino TMs by using the latest global solar neutrino 
data. We find that all elements of the TM matrix can be bounded at the same time. 
Furthermore, we show how reactor data play a complementary role to the solar 
neutrino data, and use the combination of both data sets to improve the current 
bounds. Performing a simultaneous fit of LMA-MSW oscillation parameters and TMs 
we find that 6.3 x 10^^°/iB and 2.0 x IO^^'^^ub are the 90% C.L. bounds from solar 
and combined solar -|- reactor data, respectively. Finally, we perform a simulation 
of the upcoming Borexino experiment and show that it will improve the bounds 
from today's data by roughly one order of magnitude. 
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1 Introduction 



Present solar 0, ||, ^ ^, ||, ^ and atmospheric neutrino data 0, || provide the first robust 
evidence for neutrino fiavour conversion and, consequently, the first solid indication for 
new physics. Neutrino oscillations constitute the most popular explanation for the data 
(for recent analyses see Refs. |^, |11|) and are a natural outcome of gauge theories 
of neutrino mass [0, 0. Non-zero neutrino masses manifest themselves through non- 
standard neutrino electromagnetic properties. If the lepton sector in the Standard Model 
(SM) were minimally extended in analogy with the quark sector, neutrinos get Dirac 
masses (m^) and their magnetic moments (MMs) are tiny |T^ , 



3 X lO-'V. (^) 



fll 



rzll for the electron neutrino, 

see 



where hb is the Bohr magneton. Laboratory experiments give 90% CL. bounds on the 
neutrino MMs of 1.8 x IQ-^Vb E HI and 1.3 x IO^^Vb _ 
6.8 X 10~^^^B for the muon neutrino [|T^] and 3.9 x lO^^/i^ for the tau neutrino 
also Ref. |^). On the other hand, astrophysics and cosmology provide limits of the order 
of 10~^^ to 10^^^ Bohr magnetons [^. Improved sensitivity for the electron neutrino MM 
from reactor neutrino experiments is expected, while a tritium source experiment 
aims to reach 3 x 10~^'^^b- 

It has for a long time been noticed, on quite general "naturality" grounds, that Majo- 
rana neutrinos constitute the typical outcome of gauge theories fl^ . On the other hand, 
precisely such neutrinos also emerge in specific classes of unified theories, in particular, in 
those employing the seesaw mechanism |T3[. If neutrinos are indeed Majorana particles 



the structure of their electromagnetic properties differs crucially from that of Dirac neu- 
trinos ||2^, being characterized by a 3 x 3 complex antisymmetric matrix A, the so-called 
Majorana transition moment (TM) matrix. It contains MMs as well as electric dipole mo- 
ments of the neutrinos. The existence of any electromagnetic neutrino moment well above 
the expectation in Eq. ([l|) would signal that some special mechanism — which goes beyond 
the SM — is at work. Therefore, neutrino electromagnetic properties are interesting probes 
of new physics. 

As noted in Ref. p3| the existence of Majorana transition moments leads to the phe- 
nomenon of spin-fiavour precession, when they move through a magnetic field, as might 
happen in the Sun. This was subsequently shown to affect the propagation of solar neu- 
trinos in an important way, due to the effects of matter and this has recently been 
shown to provide excellent fits to current data P3 . 



In the present paper we will take a different attitude. Motivated by the now ro- 
bust status of the LMA-MSW solution of the solar neutrino problem indicated by recent 
data |1^, and also by the prospects of its possible confirmation by the KamLAND ex- 
periment p6|, we take up this solution as the basis of our investigation. It has been shown 
in Ref. that in this case the effect of a solar magnetic field can be safely neglected. 
Therefore, our present approach is justified and complementary to that in Ref. |2^: there 
we were concerned with spin-fiavour-conversions of solar neutrinos, sensitive to the prod- 
uct of A times a solar magnetic field, while here we are concerned with the effect of TMs 
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in the neutrino detection process. Sizable TMs would contribute to the elastic neutrino- 
electron scattering in the Super-Kamiokande (SK) experiment 0], and hence, data from 
this experiment can be used to constrain electromagnetic neutrino properties 

Here we will tacitly employ the results of the latest analysis in Ref. as well as the fit 
strategies used there, to perform a fit to the global solar neutrino data in terms of LMA- 
MSW oscillation parameters and neutrino TMs, in order to derive bounds on the latter. It 
is important to clarify for which combinations of the three complex moments the bounds 
are valid. We will show that all three moments can be bounded, however, not all at the 



same level. In addition, we will show that data from reactor neutrino experiments [0, 
161 , [TtH which use elastic neutrino-electron scattering for neutrino detection are to some 



extent complementary to the solar neutrino data, so that by combining both data sets 
one can improve the bounds. We also explore the potential of the upcoming Borexino 
experiment p9|, ^ in probing the TMs, showing that it is substantially more sensitive to 
the TMs than present experiments. 

The paper is organized as follows. In Section |^ we give our theoretical framework and 
we discuss in detail which elements of the TM matrix are constrained. In Section |^ we 
present the statistical method we apply to calculate bounds on the TMs, and in Section ^ 
we describe our analysis of solar and reactor neutrino data. The resulting bounds are dis- 
cussed in Section ^ and the sensitivity of the upcoming Borexino experiment is considered 
in Section |. A brief discussion concludes the paper in Section ^ 



2 The theoretical framework 



In this paper, we adopt the following basic working assumptions: 

• Only three light neutrinos exist, doublets under SU{2) ® U{1) (no light sterile neu- 
trinos); this is well- motivated by recent global fits of all available neutrino oscillation 



data, including also those that follow from short baseline searches . 

The solution of the solar neutrino puzzle is given by the LMA-MSW solution, as 



indicated by recent global fits of solar data |T0|, [111 . 

Neutrinos have Majorana nature as expected from theory |T3|, ^| 



• Neutrinos are endowed with TMs arising from some unspecified electroweak gauge 
theory. 

Although the latter is probably the least motivated of our assumptions, there are some 
models of the electroweak interaction where the possibility of enhanced TMs has been 



discussed [13 



In experiments where the neutrino detection reaction is elastic neutrino-electron scat- 
tering, like in SK, Borexino and some reactor experiments [^5], |16|, [l^, the electromagnetic 
cross section is 

(icTpm a^vT /I 1 



dT 



(2) 
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where the effective MM square is given by |33| 



/igfj = aLA^Aa_ + a\_\\^a^ . (3) 

The electromagnetic cross section adds to the weak cross section and allows to extract 
information on the TM matrix A. In this cross section, T denotes the kinetic energy of the 
recoil electron, E^y is the energy of the incoming neutrino and the 3-vectors a_ and a+ 
denote the neutrino amplitudes for negative and positive helicities, respectively, at the 
detector. The square of the effective MM given in Eq. (H) is independent of the basis 



chosen for the neutrino state |^ . In what follows we consider both the flavour basis and 
the mass eigenstate basis. We will use the convention that a=p and A = (Aq/j) denote the 
quantities in the flavour basis, whereas in the mass basis we will use a=p and A = [Xjkj- 
The two bases are connected via the neutrino mixing matrix U: 

a_ = U^a^ , a+ = U^a+ , A = U^XU . (4) 

One usually decomposes the transition moment matrix as 

X = fi — id , (5) 

where and d are hermitian matrices. On general grounds these matrices are, in addition, 
antisymmetric and imaginary for Majorana neutrinos (see also Ref. and references 



cited therein). It is, furthermore, useful to define vectors A = (Aq,) and A = (Aj) in the 
flavour and mass basis, respectively, by 

Xap = £ai3^-Ay and Xjk = £jki-Ai . (6) 

Thus, in the flavour basis we have Ag^ = A,-, A^^ = Ag and A^e = A^. Note also that 

|A|2 = iTr(AU) ^ |A| = |A|. (7) 

This means if we are able to find a bound on |A|, we have not only constrained the TMs 
in the flavour basis but also in the mass basis. In the numerical section of the paper, this 
will be exactly our strategy. 

Let us now discuss the detailed form that the effective MM square ii^g in Eq. takes 
when we assume the LMA-MSW solution of the solar neutrino problem, denoted by /ilma- 
Numerical integrations of the neutrino evolution equations performed in Ref. showed 
that for an effective MM of 10~^^ fiB a solar magnetic field of the order of 80 KGauss has 
practically no effect on the LMA-MSW solution. The results of Ref. pS) imply that it is 



safe to neglect possible solar magnetic field effects on the evolution of the neutrino state in 
the sun in case of LMA-MSW. This follows essentially from the fact that the terms in the 
evolution Hamiltonian related to the LMA-MSW mass-squared difference and the matter 
potential are much bigger than the electromagnetic terms. Hence, helicity is conserved 
in solar neutrino propagation, so that we can set = 0, thereby eliminating the second 
term in Eq. (Rl) . Another good approximation used in fitting solar data in a three- neutrino 
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scenario [0 is to set Am^^^ ^ -^"^sol = Am^, i.e., oscillations with Am^^^ are averaged 
out. Using the parameterization @ we obtain from Eq. @ 

^^lMA = |AP - Pi!r|A3p - E A*A, . (8) 

j,k=i 

The brackets (...) in the last term in Eq. denote the average over the production point, 
earth-sun distance and zenith angle, and we have made use of the relation 

= (J = 1,2,3) (9) 

for the probability that the neutrino produced in the core of the sun as a arrives at 
the detector as a mass eigenstate uj. From Eq. (H) one learns that besides the absolute 
values of the TMs |Ai|, IA2I, IA3I only one complex phase 6 = arg(A^A2) enters into /iLMA- 
Eq. (ID can be further simplified by making use of the large mass-squared difference 
Am^ ~ 5 X 10^^ eV^ implied by the LMA-MSW solution. First, because of vacuum 
oscillations on the way from the sun to the earth the neutrino state arriving at the earth 
can be considered as an incoherent mixture of mass eigenstates Note also that 

for the LMA-MSW solution, earth matter effects are very small so that we neglect 
them in our expression for /iLMA-Q Then we find that only the diagonal elements of the 
matrix (A^A^ contribute, and as a result only the probabilities (P) appear in determining 
the detected neutrino signal. Setting S13 = sin6'i3, C13 = cos6'i3 and Ues = S13 one has, to 



a very good approximation |TT 



Pjr = 4 and PeT = C?3PeT(j = l>2) with P^r + Pe2=l- (10) 

The probabilities P^j' {j = 1, 2) are the effective 2-neutrino probabilities for the solar 
neutrino problem, where all the averages mentioned above have been taken into account. 
Using now the probabilities in Eq. ( pUj ) and the approximations mentioned above, Eq. (^) 
becomes 

/iLMA = iAir + iA2r + 4iA3r - 4 e nTiA.r • (n) 

i=i,2 



It is well known from reactor neutrino data |^ that the mixing angle Oi^ is rather 
small. Fits to the CHOOZ data and atmospheric neutrino data [|ll| show that at 3a one 
has s^3 ^ 0.05. This allows us to replace cf^ in Eq. (|ll]) by 1. Then we arrive at the final 
formula 

/^LMA = |AP - IA2P + (IA2P - |Aip) , (12) 

which will be used in the numerical section evaluating solar neutrino data. The probability 
P^i is a function of the ratio Am? / and the solar mixing angle 9 = 612. In the following 
we will drop the super-script 2v and Pel or P^e refers always to a 2-neutrino probability. 
Eq. ([T^) naturally makes no distinction between MMs and electric dipole moments. Since 

^Note, however, that the neglection of earth matter effects is not a very good approximation when 
performing a combined fit of TMs and osciUation parameters. Therefore, in this case we do take into 
account earth matter effects in the minimization of the with respect to the solar oscillation parameters. 
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we aim at constraining |A|, all elements of the TM matrix will be bounded at the same 
time. 

Now we move to the case of reactor neutrinos. There we have a pure v^. source. 
Therefore, we have a_ = and a+ = (1,0,0)"^ since in these experiments the baseline 
is much too short for any oscillations to develop. The resulting /ig^ relevant in reactor 
experiments is given as 

/z^ = |Aj2 + |A^|2. (13) 

From this relation it is clear that reactor data on its own cannot constrain all TMs con- 
tained in A, since Ag does not enter in Eq. ([131) . In order to combine reactor and solar 
data it is useful to rewrite ^ in terms of the mass basis quantities. With Eq. we 
readily derive 

Then with Eq. we obtain 

= |A|2-cVip-s^|A2|^-2sc|Ai||A2|cos(5, (15) 

where c = cos0 and s = sinO, 9 being the solar mixing angle. Further, we notice that 
the relative phase 5 = arg(A*A2) between Ai and A2 appears in addition to |A|, |Ai| and 
IA2I. 



3 Statistical method and qualitative discussion 

In the following we will use data from solar and reactor neutrino experiments to constrain 
neutrino TMs. The x^-function obtained from the data depends on the solar oscillation 
parameters tan^ 6 and Am^ as well as on the elements of the TM matrix A. Regarding the 
dependence on the oscillation parameters we will take two different attitudes. One is to 
assume that tan^ 6 and Am^ will be determined with good accuracy at the KamLAND 
experiment, and hence, we will consider the at fixed points in the tan^ 6 — Am^ plane 
(method I). In the second approach we will derive a bound on the TMs by minimizing the 

with respect to tan^ 9 and Am^ (method II). This second procedure takes into account 
the present uncertainty of our knowledge of the oscillation parameters. 

Let us describe in detail how we calculate a bound on the TMs. Our aim is to con- 
strain |A|, therefore it is convenient to consider the ^is a function of the independent 
parameters |A|, |Ai|, IA2I and 6. As discussed in the previous section, 6 appears only if 
reactor data are included. The x^-functions which we are using for the individual data 
sets (solar rates, SK recoil electron spectrum, reactor data) will be described in detail 
in the following sections. When performing the fit to the data, we find that in general 
the minimum of the occurs close or outside the physical boundary of the parameters 
|A|, |Ai|, IA2I. To take this into account we apply Bayesian techniques to calculate an 
upper bound on |A|. We minimize the with respect to |Ai|, IA2I and 6 for each value 
of |A|, taking into account the allowed region < |Aip + IA2P < |Ap: 

x2(|A|) = Min[x'(|A|,|Ai|,|A2U)l ■ (16) 
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In method I we do this for fixed values of the oscillation parameters, scanning over the 
LMA-MSW region, whereas in method II we minimize also with respect to tan^ 6 and 
Am^ in Eq. ([T6|). Then the is transformed into a likelihood function via 

C oc exp i-^pd) ■ (17) 

Now we use Bayes' theorem and a fiat prior distribution p(|A|) in the physically allowed 
region, p(|A|) = G(|A|), to obtain a probability distribution for |A|: 

/:(|A|)e(|A|)c^|A| 

S^C{\K\')d\K\' ■ ^ ^ 

An upper bound ha on |A| at a C.L. a is given by the equation 



/(|A|)rf|A| 



/(|A|)rf|A| = a. (19) 

Let us consider in more detail the minimization with respect to |Ai| and IA2I in Eq. ([T6|) 
in the case of solar data without reactor. As we will show later there is no evidence for a 
non-zero /Ilma data. Therefore the minimum of the for a given | A| will occur if 

/^LMA is minimal. Departing from Eq. (|T^ it is easy to show that 



2 / |A|2p,i for Pel < 0.5, 

Mm[/iLMA]= ,.,2/1 D ^ D . n. (20) 



The minimum occurs at 



|A| 



Pel > 0.5. 



IA2P for Pel < 0.5, 

|Ai|2 + |A2|2 for Pel = 0.5, (21) 
|Ai|2 for Pel > 0.5. 



From Eq. ( PU] ) follows that the bound on |A| is strongest if Pei = 0.5 because in this 
situation /Ulma is maximal. In Fig. |l] we show contours of constant Pei in the tan^ 9 — 
Am? /Ey plane. For definiteness, the probabilities in the figure are obtained by performing 
the averaging over the production distribution inside the sun for the '^Be fiux most relevant 
for Borexino. However, the probabilities for the ^B fiux relevant for SK are very similar. 

The shaded region in Fig. |I| is the region relevant for Borexino, assuming a mass- 
squared difference in the range 10~^ eV^ < Am^ < 10~^ eV^, whereas for SK the region 
below the dash-dotted line is most important, due to the higher energy of the ^B neutrinos. 
We can read off from the figure that in most part of the SK region Pei is very small. 
According to Eq. ( PT| ) this means that the sensitivity of SK for |A| is limited because of 
the small Pei. In contrast, we expect a much better sensitivity of Borexino, because, in 
a large part of the relevant parameter space in this case, Pei is close to the optimal value 
of 0.5. 

The arguments presented here are valid as long as the depends only on the pa- 
rameter combination given in Eq. (|12D- The data from reactor experiments depend on a 
different combination of TM parameters (see Eq. (|15[)) implying that the combined analysis 
of SK and reactor data will potentially lead to more stringent limits than the SK data 
alone. 
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tan e 

Figure 1: Contours of constant Pel- The shaded region is relevant for Borexino, whereas 
the region below the dash-dotted line is relevant for SK. 

4 Analysis of solar and reactor neutrino data 

In this section we describe briefly our analysis of solar neutrino data and the data which we 
are using from reactor neutrino experiments. Neglecting the electromagnetic contribution 
to the small elastic scattering signal in SNO, one notices that the neutrino TMs contribute 
directly to the observed signal only in SK. However, the uncertainty in the flux leads to 
correlations between SK and all other experiments. Therefore, even for fixed oscillation 
parameters (method 1), also other experiments give some information on TMs and it is 
important to include them in the analysis p8| . 

We divide the solar neutrino data into the total rates observed in all experiments and 
the SK recoil electron spectrum (with free overall normalization). Before combining rates 
with SK spectra, we derive bounds on | A| for each of these cases separately. Subsequently 
we consider the reactor data by themselves and in combination with the global sample of 
solar data. 
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CI [0 


2.56 ±0.16 ±0.16 SNU 


Sage [0 


70.8 to SNU 


GaGNO § 


73.3 ±4.7 ±4.0 SNU 


SKi 


0.465 ± 0.005 t{];g|^ 


SNO CC [| 


0.3485 ±0.012 ±0.018 


SNO NC § 


1.008 ±0.087 ±0.097 



Table 1: Solar neutrino rates used in our analysis. For SK and SNO the given number is 
the ratio of the observed flux and the flux predicted by the BPOO SSM . 



4.1 Solar rates 

In this analysis we use the event rates measured in the Chlorine experiment Homestake 
lU (CI), the Gallium experiments Sage 0, Gallex and GNO 0] (GaGNO), the total event 
rate of SK |Q based on the 1496 days data sample and the recent result of the Sudbury 
neutrino observatory for the charged current (SNO CC) and neutral current (SNO 
NC) solar neutrino rates. For the analysis we use the x^-function 



[22) 



IS 



Here the indices j, k run over the 6 solar neutrino rates (CI, Sage, GaGNO, SK, SNO CC 
SNO NC), Dj are the experimental rates, Rj are the theoretical predictions and VJ^^"^ 
the covariance matrix. 

The data we are using are given in Table 0. The rates of SK and SNO are the ratio of 
the measured flux and the flux predicted by the standard solar model (SSM) BPOO 



and no disappearance of solar u^, ($g = 5.05 x 10^ cm ^s ^). Note that the SNO CC 
and NC fluxes as given in Ref. i have been obtained from a flt to the SNO spectral 
data assuming no distortion of the neutrino energy spectrum, i.e., a survival probability 
constant in energy, which in general is not realized in the case of neutrino oscillations. 
However, in the case of the LMA-MSW solution it is a very good approximation to use 
the fluxes given in Ref. 0. 

The theoretical predictions Rj which appear in Eq. (0) are calculated as described 
in Ref. However, the rate -Rsk of the SK experiment includes an extra contribution 
from the electromagnetic scattering. This rate is given by 



-Rsk 



/o» dE, ME.) {pUEu)ctY{E.) ± [1 - Pee{E.)K{E.) + f^luAiEuXJE,)} 



J,°°dE,MEu)cT^{E, 



(23) 



where (j)^{Ey) is the shape of the B flux normalized to 1, which we take from Ref. |3 
and 



(y{E, 



^-^^ jj^ dT R{T\T) with a = ar,a;r,^L 



(24) 
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We use the SM weak cross sections cr^ including radiative corrections , /Ulma given 
in Eq. (|l|), and 

dV mlfil \T' Ej ■ ^ ' 

Here E,y is the neutrino energy, while T' and T denote the true and the measured recoil 
electron kinetic energies, respectively. The integration in T is from the SK threshold 5 
MeV (total electron energy) up to 20 MeV. The maximum kinetic energy of the recoil 
electrons is given by Tj^axiE^) = 2E'l/{2Ei, + m^). The resolution function R(T',T) is 
taken from Ref. ^ . 



The covariance matrix VJ^^'^ in the of Eq. (^) takes into account the experimental 
errors and theoretical uncertainties from detection cross sections V^*^^ and SSM predictions 
Vp^^ added in quadrature |]4T[| , 



Vjr = + V,-f + • (26) 

The experimental errors for different experiments are uncorrelated. For j, k = CI, Sage, 
GaGNO, SK we have V^^^ = Sjk<7j, where the are the experimental uncertainties, 
calculated by adding in quadrature the statistical and systematic uncertainties given in 
Table |I]. However, the statistical and systematic errors for the SNO CC and NC rates are 
strongly correlated, with p*^**^* = —0.518 (the same as the correlation for the CC and NC 
day-night asymmetry, given in Ref. 0) and p^^^ = —0.508 (which we derive from Table II 
of Ref. 0). The second and third terms in Eq. (^) take into account the uncertainties in 
the detection cross sections and in the SSM predictions of the neutrino fluxes, respectively. 



For details see Refs. H, HTH and references therein. 



4.2 The Super-Kamiokande recoil electron spectrum 

In this section we consider the shape of the SK recoil electron spectrum. The electromag- 
netic contribution to the neutrino-electron scattering cross section leads to a substantially 
different spectrum of the scattered electrons than expected from the SM weak interaction. 
Therefore, the SK spectral data provide a useful tool to constrain neutrino TMs . 



We perform a fit to the latest 1496 days SK data presented in Ref. [Q. There the 
event rates are given in 8 energy bins ranging from 5 to 20 MeV. The energy bins 2 to 7 
are further divided into 7 bins of zenith angle, which makes up a total of 44 data points 
Di. We define the x^-function 

44 

= E («^^ - D,){Vr'T\aR, - D,) . (27) 

The theoretical prediction i?j(tan^ 6*, Am^, Plma) i-th bin is calculated as in 

Eq. (p3D with the integration interval [Ti , T2] in Eq. (p^) chosen according to the energy 
interval of the given bin. The covariance matrix V^J'^'^^ in Eq. (^) contains statistical and 
systematic experimental errors. These are obtained from Ref. taking into account the 
correct correlations of systematic errors. The of Eq. (^) is minimized with respect to 
the normalization factor a in order to isolate only the shape of the spectrum. 
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Figure 2: Ax^ of the reactor experiments Rovno |T5[ and Bugey [|T^ and the total Ax^gj^j,^Qj, 
function of the effective MM. 



4.3 Reactor data 



Data from neutrino electron scattering at nuclear reactor experiments ||T5| , p!q , p!7| can 
be used to constrain the combination of TMs given in Eq. (pISl). Here we use data from 
the Rovno nuclear power plant |T^ and from the Bugey reactor |T^. To include this 
information in our analysis we make the following ansatz for the x^-function: 



X: 



reactor 



obs 



:28i 



The sum is over the two experiments Rovno and Bugey, A^*j,g is the observed number of 
events with the one standard deviation error a\ A^^eak is the number of events expected 
in the case of no neutrino TMs (only the standard weak interaction) and N^^{fiB.) is the 
number of events due to the electromagnetic scattering of neutrinos with an effective MM 
/iR. We can write the latter as 



10-1 Vs. 



(29) 



T5[ the values ^"0^°/"° = 41, 



For the Rovno reactor we obtain from Table I of Ref. 
^Rovno _ 2Q QT^(\ N^°^^" = 32. Furthermore, in this table the number of events due to the 
electromagnetic scattering of neutrinos with an effective MM /xr — ^ in-io, 



iV^^'^°(2 X 10 I'^/ifi) = 54. Using this last information the constant in Eq. 



2 X 10 fiB is given: 
can be 
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Figure 3: Contours of the 90% C.L. bound on |A| in units of 10~^^^b from solar data. 
The shaded region is the 3a LMA-MSW region from Ref . , the best fit point is marked 
with a star. The crosses labeled A, B, C, D are the test points considered in Table ^ 



determined as = 54/4. Results of the MUNU experiment performed at the Bugey 

reactor have been presented recently at the Neutrino 2002 conference fl^. These lead to 
the bound /iR < 1.3 x 10~^°/iB at 90 % C.L., if an energy cut of Ee > I MeV is used. From 
Ref. [0 we get (in events per day) N^^f^ = 0.37, a^'^s'^y = 0.25 and A^^^af ^ = 0.45. The 
constant in Eq. (^) can be calculated by making use of the 90% C.L. bound cited above: 
we solve the equation xLgey(/^R = 1-3 x IO-^Vb) = 2.7. This leads to C^^^s^^y = 0.196. 
In Fig. D we show the Ax^ of the two experiments and also the sum. Note that with 
our ansatz Eqs. (pSD and (|^) we can reproduce to good accuracy the 1 a upper bound 
/xr < 1.5 X 10^^°/x^ from the Rovno data [|15| . 
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5 Bounds on |A| from solar and reactor data 

First we discuss the bounds on | A| from solar data alone. Fixing the oscillation parameters 
at the current best fit point tan^ 6 = 0.44, Am^ = 6.6 x 10~^ eV^ [^] we obtain with the 
Bayesian methods described in Section ^ the 90% C.L. bound 

|A| < 4.0 X lO'^^fiB (best fit point, solar data). (30) 

However, such a bound substantially depends on the values of the neutrino oscillation 
parameters. In Fig. ^ we show contours of the 90% C.L. bound on |A| in the tan^ 6 — Am'^ 
plane. We find that the bound gets weaker for small values of Am^, whereas in the upper 
left part of the LMA-MSW region a bound of the order 2 x 10~^^^b is obtained. 

By combining solar and reactor data we obtain considerably stronger bounds. At the 
best fit point we get at 90% C.L. 

|A| < 2.0 X 10^^^ fiB (best fit point, solar + reactor data). (31) 

In Fig. ^ we show the contours of the bound in the tan^ 6 — Am^ plane for the combination 
of solar and reactor data. By comparing Fig. § with Fig. § we find that reactor data 
drastically improve the bound for low Am^ values. 

In order to better understand these results we give in Tab. | bounds on |A| for five 
characteristic points in the tan^ 9 — /S.m? plane, including the current best fit point, which 
is labeled with a star. The other points are marked in Figs. ^ and ^ by crosses and are 
labeled A, B, C, D, and the corresponding values of the oscillation parameters are given in 
Tab. ^. Besides the bound on |A| we give also bounds on the individual |Aj| {i = 1, 2, 3). 
These bounds are obtained by setting the other two Aj to zero. We show the bounds for 
the different data samples: solar rates, the shape of the SK spectrum with free overall 
normalization (SK spect), solar rates combined with the SK spectrum (solar global), 
reactor data only, solar data combined with reactor data (sol + react). 

Notice that, in contrast with the bound on |A| which, thanks to Eq. (|^), is valid in 
any basis, the bounds on individual |Aj| do not have a basis-independent meaning. The 



bounds on |Aj| given in Tab. |^ refer to the mass-eigenstate basis, also used in Ref. [23 . 

By comparing the first and the second column of Tab. ^ we notice that both solar 
rates and SK spectrum are important: although the bound from the SK spectrum is 
slightly stronger, the one which follows from solar rates is in general of the same order. 
The bound we obtain by combining rates and SK spectrum is displayed in the third 
column. Concerning the reactor data, clearly these bounds do not depend on Am^, and 
the bound on IA3I is even completely independent of the neutrino oscillation parameters 
(see Eq. (|^)). As mentioned above, reactor data on their own cannot be used to constrain 
|A|. 

In the upper parts of the LMA-MSW region, the solar data alone give already a strong 
bound on |A|. From the bounds given in Tab. ^ at point A we see that for such high 
values of Am^ all |Aj| are strongly constrained by solar data. This is due to the fact that 
there the probability Pel relevant in SK is close to the optimal value of 0.5 (see Fig. |1|). 
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Figure 4: Contours of the 90% CL. bound on |A| in units of 10~^°/xb from combined 
solar and reactor data. The shaded region is the 3a LMA-MSW region from Ref. 0, the 
best fit point is marked with a star. The crosses labeled A, B, C, D are the test points 
considered in Table |^. 

In this parameter region the combination with reactor data does not improve the bound 
significantly. This is also true for the region of small tan^ 6 around point C. 

The behavior of the bounds for low Am^ values can be understood by looking at point 
B. In this Am^ region the probability Pel is very small for the neutrino energies relevant 
for SK, as seen in Fig. |l[ Therefore, solar data give only a very weak constraint on IA2I as 
can be seen from Tab. ^ and the bound on |A| is dominated by this weak bound on IA2I 
(see Eq. (|21|)) leading to the rather poor bound of 11.6 x 10^^°//^ from solar data. In this 
case the combination with reactor data improves the bound drastically to 2.3 x lO^^^fiB- 

Up to now we have calculated bounds on neutrino TMs for fixed values of the oscillation 
parameters tan^ 6' and Am^ (method 1 described in Section |^). These results will be 
especially useful after KamLAND will have determined the oscillation parameters with 
good accuracy. In the following we change our strategy and minimize the for each 
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Table 2: 90% C.L. bounds on |Aj| {i = 1,2,3) and on |A| in units of 10"^°/!^ at fixed 
points in the tan^ 9 — Am^ plane from different data samples. The hyphen indicates that 
lAI cannot be bounded from reactor data alone. 



value of |A| with respect to tan^ 6' and Am^ (method II). This will lead to a bound on 
|A| taking into account the current knowledge concerning the oscillation parameters. To 
this end we make use of the global solar neutrino data (including earth matter effects) 
as described in Ref. Q in order to obtain the correct behavior as a function of tan^ 9 
and Aim?. Only in the expression for /Xlma neglect earth matter effects as described 
in Section |^, since they are known to be small in the LMA-MSW region |3^. Performing 
this analysis we obtain the following bounds at 90% C.L.: 



^ I 6.3 X 10 ^^jJiB (unconstrained, solar data) 

2.0 X IQ^^^jiB (unconstrained, solar + reactor data) 
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Together with Figs. |^ and ^ the bounds given in Eq. (^) constitute an important 
result of this paper. They show that, assuming the solar LMA-MSW solution, current 
solar neutrino data can be used to constrain all elements of the Majorana neutrino TM 
matrix. A combination with data from reactor experiments significantly strengthens the 
bound on lAI. 



6 Simulation of the Borexino experiment 



Here we investigate the sensitivity of the Borexino experiment p9|, |3^ to neutrino TMs. 
This experiment is mainly sensitive to the solar ^Be neutrino flux, which will be measured 
by elastic neutrino-electron scattering. Therefore, Borexino is similar to SK, the main 
difference is the mono-energetic line of the ^Be neutrinos, with an energy of 0.862 MeV, 
which is roughly one order of magnitude smaller than the energies of the ^B neutrino flux 
relevant in SK. 



6.1 The Borexino x^-function 

To estimate the sensitivity of Borexino we consider the following x^-function: 

xLexino = E (^f - mV!;-T\Nf^ - D,) . (33) 

Here A^**^ is the theoretical prediction for the number of events in the electron recoil 
energy bin i, and Di is the (hypothetical) observed number of events. Borexino will 
observe recoil electrons with kinetic energy in the range 0.25 to 0.8 MeV. As a realistic 
example we adopt A^bins = 8 bins in electron recoil energy and an energy resolution of 



0.058 MeV |g9|, |30|, We have checked that our results are essentially independent of 
the exact value of the energy resolution. 

The theoretical prediction for the number of events in the bin i after Ayears years of 
Borexino data taking is given by 

with X = th, weak, em, bg, and fij'^^^^'^'^^^s jg f_]^Q number of events per day from weak 
scattering, electromagnetic scattering and the background, respectively. If not stated 
otherwise we consider a running time of Ayears = 3 years. In order to estimate the 
sensitivity of Borexino for neutrino TMs we assume that the data are generated by neutrinos 
without TMs: 

= NJ""^ + N^^ . (35) 

Hence, we obtain 

xLrexino = E (^J"^^^)-' . (36) 

The minimum of this occurs for /xlma = and is always zero. With the of Eq. (|36| ) 
we calculate a bound on |A| at a given CL. This bound corresponds to the maximum 
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flux 


'Be 


pep 


i3]\j- 15q background 


events/ day 


43.3 


2.0 


4.0 5.5 19 



Table 3: Expected number of events per day in Borexino for the kinetic energy of the 
recoil electrons in the range 0.25 to 0.8 MeV resulting from different components of the 
SSM solar neutrino flux without neutrino conversion (taken from Table 4 of Ref. [pOp . 
Also shown is the expected number of background events per day. 



allowed value of |A| which cannot be distinguished from |A| = 0, and is therefore a 
measure for the obtainable sensitivity to |A| at Borexino. 

The numbers of events due to the weak and electromagnetic scattering are calculated 

by 



^weak 



f ''^ 

/•CO 

-^E"^// C^E.0y/i^MA[<m]., (37) 
/ 

respectively. The effective cross sections in the bin i are calculated like in Eqs. (|2^) and 
(^), by choosing the lower and upper bounds of the recoil energy interval Ti and T2 
according to the given bin i. We sum over the four most important neutrino fluxes (/ = 
Be, pep, N, O) in Borexino (see Ref. ||3^). For the mono-energetic lines we have (pBciEu) = 
6{E^ — 0.862 MeV) and 4>pcp{E^) = 5{Ey — 1.442 MeV), and we use the normalized spectra 
4>-^{Ei^) and (f)o{E^) given in Ref. [^. The absolute values of the neutrino fluxes $° from 
Ref. 1^ are given by = 0.429, <l>°,p = 0.014, $^ = 0.055, = 0.048 in units of 10^° 
cm~^s~^. Note that from the two ''Be lines the line at 0.38 MeV, which constitutes 10% 
of the total ^Be flux, does not contribute to the signal in Borexino. Therefore the value 
for given above has been obtained by multiplying the value given in Ref. with 
0.9. 

In Table |] we show the expected number of events per day in Borexino for the kinetic 
energy of the recoil electrons in the range 0.25 to 0.8 MeV, resulting from different solar 
neutrino fluxes, for the standard solar model without neutrino oscillations. We do not 
consider other neutrino fluxes, because they contribute to the signal with less than 0.2 
events per day. We flx the normalization constant A/" in Eq. ( |37D in such a way that the 
number of events per day from the ''Be line for Pee = 1 and /zlma = is 43.3. Then we 
can reproduce the number of events arising from the other fluxes as given in Table ^ to 
good accuracy. 

The number of background events n^^ is obtained in the following way. We read off 
the shape of the background as a function of the electron energy from Fig. 15 of Ref. 



and normalize this function to 19 background events per day in the energy range 0.25 to 
0.8 MeV p9|. Then we can calculate n^^ by integrating over the corresponding energy 



interval for each bin. 
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tan e 

Figure 5: Contours of the 90% C.L. bound on | A| after 3 years of Borexino data-taking in 
units of 10~^^ fiB- The current best fit point is shown by the star, and the shaded region 
is the allowed LMA-MSW region at 3a, from Ref. The dotted line corresponds to 
Pel = 0.5 for ^Be neutrinos. 

We use the following covariance matrix in Eq. (^): 

V^T"" = + E N^f^N.j, f: af,,af,, (AlnX^)^ + iV^iV]'^(A In N^^^f . (38) 

/l,/2 13=1 

The first term is the statistical uncertainty, taken as the square-root of the predicted 
number of events. The second term describes the uncertainty in the solar neutrino fluxes. 
Here Nij are the contributions of the individual solar neutrino fluxes to the total event 
numbers in each bin: A^**^ = J2f ^^if + ^i^^? and the sum is over the four fluxes relevant in 
Borexino (for definition and values of the quantities a//? and see Ref. The last 

term takes into account a systematic uncertainty in the number of background events. 
We assume full correlation between the bins and that the relative uncertainty is the same 
for all bins: AN^^/N^'^ = A\nN^^. We adopt a standard value for our calculations of 
AlniV^g = 10%. 
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Am^ = 6.6x 10 ^eV^ 




K. =8 

bins 



N 



Figure 6: The left panel shows the 90% CL. bound on | A| after 3 years of Borexino data- 
taking as a function of the number of recoil electron energy bins for two values of the 
error on the background. The right panel gives the 90% CL. bound on |A| as a function 
of Borexino running time for AlniV^^ = 10%. 



6.2 Sensitivity on |A| of Borexino 



Using the statistical method described in Sections ^ and [GT^ we obtain for the current 
best fit point Am^ = 6.6 x 10~^ eV^, tan^ 9 = 0.44 the upper bound (sensitivity) 



A| < 0.28 X IO^^Vb at 90% CL. 



(39) 



after three years of Borexino data taking. This bound is about one order of magnitude 
stronger than the bound from existing data. We have checked that a combined analysis 
of Borexino with existing data (solar as well as reactor data) does not improve the bound 
of Eq. (^). The bound depends on the actual value of the oscillation parameters; in 
Fig. § we show contours of the 90% CL. bound in the tan^ 6 — Am^ plane. The strongest 
attainable limit is roughly 0.24 x IO^^^/^b- In agreement with the discussion we gave in 
connection with Eq. (pT]) we find that the strongest bound is obtained when P^i = 0.5, as 
shown by the dotted line in Fig. ^ 

We have performed several tests concerning our assumptions for Borexino. Fig. ^ (left 
panel) shows the dependence of the bound on the number of bins in the recoil electron 
energy. We can see from this figure that a measurement of the total rate alone (A^bins = 1) 
gives a rather weak bound. However a few bins are already enough for a considerable 
improvement of the bound on |A|. It is interesting to see that the error on the background 
is only important for few bins, and becomes completely unimportant for more than 3 or 
4 bins, because the main information comes from the spectral shape. 

In the right panel in Fig. ^ we show the bound as a function of the Borexino running 
time. Here we adopt again our nominal value of 8 energy bins. After one year a sensitivity 
of 0.4 X 10~^°/is can be obtained with Borexino, which is already about one order of 
magnitude stronger than the current bound. After three years of data-taking only a 
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minor improvement of the bound is possible, mainly due to the uncertainty in the solar 
neutrino fluxes.0 



7 Conclusions 

In this paper we have presented stringent bounds on electromagnetic Majorana transition 
moments (TMs). Such TMs — if present — would contribute to the elastic neutrino-electron 
scattering signal in solar neutrino experiments like Super-Kamiokande or the upcoming 
Borexino experiment, as well as in reactor neutrino experiments. Motivated by the robust 
status of the LMA-MSW solution of the solar neutrino problem indicated by recent data, 
and also by the prospects of its possible confirmation by the KamLAND experiment, we 
have taken this solution as basis for our investigation. In this context we have clari- 
fied which effective magnetic moments are probed in solar and reactor experiments — see 
Eqs. ([T^) and ([T3|), respectively. Using most recent global solar neutrino data we have 
performed a fit in terms of the oscillation parameters and the elements of the complex TM 
matrix A of three active Majorana neutrinos. Taking into account the antisymmetry of 
the TM matrix by the parameterization Aj^ = EjkiAi, we have shown that solar neutrino 
data allow to constrain |A|, i.e., all elements of the TM matrix are bounded at the same 
time. We want to stress that the bounds on |A| hold in any basis, because |A| is an 
intrinsic neutrino property; in terms of TMs it is given by 

|A|2 = |Ae;.|' + |A^,|2 + \\re\' = l^u]^ + [Assl' + M (40) 

in the flavour and mass-eigenstate bases, respectively. On the other hand, the bounds on 
the individual |Aj| {i = 1, 2, 3) in Tab. § refer only to the mass-eigenstate basis. 

A fit to the global solar neutrino data leads to the bound |A| < 6.3 x 10~^^ fiB at 
90% CL. We have also considered the role of reactor neutrino data on neutrino TMs, 
shown to be complementary to solar neutrino data. A combined fit of reactor and solar 
data leads to significantly improved bounds: at 90% CL. we get |A| < 2.0 x lO^^^fi^. 
In the very near future the KamLAND experiment will crucially test the LMA-MSW 
solution. If KamLAND confirms the LMA-MSW solution then a precise determination 
of the parameters tan^ 6 and Am^ might be possible. This motivated us to scan the 
tan^ 6 — Arn^ plane and to calculate the corresponding bound on |A| in each point. These 
results are shown in Figs. |^ and § for solar data and solar + reactor data, respectively. 

We have also investigated the potential of the upcoming neutrino-electron scattering 
solar neutrino experiment Borexino to constrain neutrino TMs. Performing a detailed 
simulation of the experiment we find that it will improve the bound on |A| by about one 
order of magnitude with respect to present bounds. 

Last, but not least, let us mention that alternative solutions to the solar neutrino 



problem based on non-standard neutrino matter interactions [E3|, which may arise in 



models of neutrino masses P3, can also be constrained in an analogous way |H2[. If, as 



seems likely, the LMA-MSW solution is finally borne out, the improved determination of 



^Note that we do not include any experimental systematic uncertainties in our Borexino analysis. 
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the oscillation parameters expected, say at KamLAND, will provide an extremely useful 
stepping stone for probing other features of neutrino physics beyond the Standard Model. 
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